Accepted for publication in the Astrophysical Journal Letters 
Preprint typeset using LfTgX style emulateapj v. 1 1/12/01 



GEMINI NEAR INFRARED SPECTROGRAPH OBSERVATIONS OF A RED STAR FORMING GALAXY 
AT Z = 2.225: EVIDENCE FOR SHOCK-IONIZATION DUE TO A GALACTIC WIND 1 

P. G. VAN DOKKUM 2 , M. KRIEK 3 , B. RODGERS 4 , M. FRANX 3 , AND P. PUXLEY 4 
Accepted for publication in the Astrophysical Journal Letters 

ABSTRACT 

Recent studies have shown that /^-luminous galaxies at 2 < z < 2.5 have high Ha line widths and high 
[Nil] A6583/Ha ratios. If these lines arise from photo-ionization by massive stars in Hll regions these results 
imply that massive, metal rich galaxies exist at high redshift. Here we investigate the ionization mechanism in 
a galaxy with K s = 19.1 and z = 2.225 in the Chandra Deep Field South, using the new Gemini Near Infrared 

' Spectrograph (GNIRS). GNIRS' cross-dispersed mode gives simultaneous access to the entire 1 - 2.5 /im 

wavelength range, allowing accurate measurements of line ratios of distant galaxies. From the ratio of Ha/H/3 
we infer that the line emitting gas is heavily obscured, with E(B — V) = 0.81q2- The reddening is higher than 
that inferred previously from the UV-optical continuum emission, consistent with findings for nearby star burst 
galaxies. We find that the galaxy has Seyfert-like line ratios, [Nil] /Ha « 0.6 and [Olll] A5007/H/3 w 6, which 
QJ ' can be caused by photo-ionization by an active galactic nucleus (AGN) or shock ionization due to a strong galactic 

wind. Although we cannot exclude the presence of an AGN, the lack of AGN spectral features in the rest- 

• frame ultraviolet, the consistency of radio, X-ray, and rest-frame UV star formation indicators, the fact that the 

[O III] /H/3 ratio remains high out to ~ lOkpc from the nucleus, and the observed gas kinematics all argue for the 
t— I . wind hypothesis. Similar shock-induced ionization is seen in nearby star burst galaxies with strong winds, such 

as NGC 1482 and NGC 3079. The evidence for shock ionization implies that measurements of metallicities and 

, dynamical masses of star forming z > 2 galaxies should be regarded with caution, especially since the existence 

of strong galactic winds in these objects is well established. Based on Sloan Digital Sky Survey data for nearby 
galaxies and the limited data available at high redshift we speculate that the effects of shocks may correlate with 
dust content. The results presented here demonstrate the importance of measuring the full rest-frame optical 

■ spectra of high redshift galaxies, and showcase the potential of GNIRS for such studies. 
Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation 

O 1. INTRODUCTION galaxies at z > 2 is well established (e.g., Franx et al. 1997; 

XT • ■ ,* TTr> s t • j , U1 • e Pettini et al. 1998, 2001). The ionization source cannot be 

O ■ Near-intrared (NIR) spectroscopy provides valuable infor- . , r ' _ , TT ' , , , , 

. .. . „ , .• j i • ^ c i • » constrained from [N II] / Ha alone (although values larger than 

5— i mation on the stellar populations and kinematics of galaxies at _ . . , L , . , b . . . b . TT 

» ■ ^ .. . f . £ .. , .... , w 0.5 point to other mechanisms than photo-ionization in H II 

c/3 z > 2, as it gives access to rest-frame optical emission lines that \ , , ,. , „, , , r , , 

—> , , . ,. , . . , . , , TT . XITri regions), and the standard method of spectral classification 

■ have been studied extensively in the local Universe. NIR spec- f. , • r ■ , „, 

• • ■ t c ii j Mir. i ^ j i • » <-> ^ „ i u relies on combinations of line ratios such as [N II] / Ha and 
tra of optically- and NIR-selected galaxies at 2 < z < 3 have A5007/H/3 (e Baldwin Philli s & Terlevich 198L 

constrained or measured their star formation rates, metallici- v ill &Ot b k 1987") ^ ' 1 ^ S ' er evlc 
ties, wind velocities, and masses (e.g., Pettini et al. 2001; Erb . , ,. , . . ,. , , , . r . 

s_j ■ t i on m n ,, t , ,„„, cu i . 1 onn/i\ Measuring these line combinations at high redshift is difri- 

H , et al. 2003; van Dokkum et al. 2004; Shapley et al. 2004). , , , ,. . , , r 

C3 . ^ ■ • • c i- . , 4, cult due to the limited wavelength coverage of most existing 

Among the most surprising results is the finding that apparently XTT _ , , , b 

t , . .5 j j ,. c °nnni -l j NIR spectrographs on large telescopes. In this Letter, we use 
mature galaxies with Ha line widths of a > 200 km s and , r , b r . . ,°_ . . . % T _ 
rxT TT n \/coi ,n t - c ^ n 1 i a ' **j * r, c / the newly commissioned Gemini Near Infrared Spectrograph 
[Nil] A6583/ Ha ratios of > 0.3 already existed at z~ 2.5 (van , , ,. . „ T ^ , , „ & . 
n n . , OAA/l cu l * i iriA^ (GNIRS) to study the line ratios of a A -selected star forming 
Dokkum et al. 2004; Shapley et al. 2004). v , ' _ JL. . , _ , T „. _ , „f 
rp, • t . .. . tu t . t tU ■ ■ galaxy at z = 2.225 in the Chandra Deep Field South (CDF-S). 
The interpretation rests on the assumption that the emission J^™,, , , . , , , . , 
,. c u » • • u • * ■ „„ GNIRS can be operated in cross-dispersed mode, which pro- 
lines arise from photo-ionization by massive stars in H II re- , r . . „ . \ , 

i T3 u * i onm cu i . 1 iaa/in t duces spectra over the entire 1 urn- 2.5 urn wavelength range 

gions (see, e.g., Erb et al. 2003; Shapley et al. 2004). In v . . . _ , , . , 6 

~, , i it • a- • ,m i j. • •. at a resolution of R ~ 1800. GNIRS thus allows simultaneous 

the local Universe this assumption is valid for the maionty „ . r „■ 

c . n ■ , • , . . n ,, rn. rMt , n ,„ y measurement of rest-frame optical lines falling in the 7, H, and 

of star forming galaxies, but not for all. The [N II] / Ha ratio „ , , , , • , • , , , , , 

u a ju - • s u a/^m A bands through identical atmosphere and aperture on the sky, 

can be influenced by photo-ionization by an AGN continuum, ,, 6 , «• . , • • 

, , • • .. ii . ^ ... c ^ a ■ and hence accurate and efficient determinations of line ratios, 

shock ionization caused by the interaction of outflowing gas 

with the ambient interstellar medium, photo-ionization by hot, 

high metallicity stars in a dense environment, and other pro- 2. SPECTROSCOPY 

cesses (see, e.g., Veilleux et al. 1994, and references therein). , ™, o , „ Tm 

The potential effects of shock ionization are particularly rele- ^ Y C ® FS ~? 95 T 2004 ^ptember 2-3 UT 

vant, as the existence of large scale outflows in star forming wlth GNI ^ S on C ^™ m "^° Uth ' , Under m ° Stl , y cle f,_ skles wlth r 

image quality w 0. 7. The galaxy was selected because of 
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Fig. 1. — Cross-dispersed low resolution spectrum of CDFS-695 at z — 2.225 (top panel), smoothed with a 200 km s _1 boxcar filter. Grey bands 
indicate regions of strong atmospheric emission or absorption. The familiar rest-frame optical emission lines are readily identified, showing the 
potential of GNIRS for measuring redshifts and line ratios in the near-IR. The bottom panels show (unsmoothed) details of the spectrum, centered 
on emission lines. 



its ^-magnitude and red J — K s color. It has K s = 19.1 and 
J — K s = 2.3 (S. Wuyts et al., in preparation), and is one of 
the most luminous members of the Distant Red Galaxy (DRG) 
class (Franx et al. 2003; van Dokkum et al. 2003) in the CDF- 
S. After measuring the redshift with GNIRS we realized that 
the galaxy is object 5 in the Daddi et al. (2004) sample of K- 
selected star forming galaxies. We note that Daddi et al. (2004) 
give a slightly bluer color of J — K s = 2.15 for this object; we 
cannot determine the source of this small discrepancy. 

GNIRS was operated in cross-dispersed mode, using the 
short camera with the 32 lines mm -1 grating and the 0"68 
slit. Between individual ~ 300 s exposures the objects was 
moved along the 6. "2 slit in an ABA'B' pattern, with A' and 
B' small (0"3) perturbations of A and B. The total exposure 
time was 5500 s. The data reduction used a combination of 
standard procedures and custom software, and is described in 
M. Kriek et al., in preparation. The reduced spectrum is shown 
in Fig. 1. The continuum is clearly detected along with famil- 
iar rest-frame optical emission lines, which provide a redshift 
Z = 2.225. This high quality spectrum demonstrates the poten- 
tial of GNIRS for studies of faint objects, and specifically for 
measuring redshifts of ^-selected galaxies. 

Line fluxes were determined from the residual emission after 
subtraction of a linear fit to the continuum, with errors esti- 
mated from repeat measurements in empty regions. Equivalent 
widths were determined from a direct comparison of the line 
fluxes to the continuum fits, and the line fluxes were calibrated 



by comparing the measured total K s magnitude to a synthetic 
magnitude determined from the GNIRS spectrum. This method 
assumes that the line emitting gas has the same distribution as 
the continuum emission, and we estimate that the systematic 
error is w 20%. Finally, line luminosities were calculated us- 
ing n m = 0.3, n A = 0.7, and H = 70kms _1 Mpc" 1 . Fluxes, 
luminosities, and rest-frame equivalent widths are listed in Ta- 
ble 1 . Errors do not include a ~ 20 % uncertainty in the overall 
normalization of the spectrum (which cancels for line ratios). 

3. REDDENING 

The reddening of the line-emitting gas can be derived from 
the Balmer decrement, after correcting for underlying contin- 
uum absorption at Ha and H/3. Stellar populations with ages 
0.5 - 2Gyr (Forster Schreiber et al. 2004; Daddi et al. 2004) 
have Balmer absorption equivalent widths of 4± 1 A. The cor- 
responding corrections to the measured line luminosities are a 
factor of 1.04 for Ha and a factor of 1.3 for H/3, and the cor- 
rected ratio Ha/H/3 = 6.6^j \. Assuming an intrinsic ratio of 
2.76, which is valid for Hll regions and also holds to good 
approximation for shock ionized regions (Binette, Dopita, & 
Tuohy 1985), the reddening E[B - V) = 0.8+^2. This large 
value is consistent with previous indirect determinations for 
DRGs (van Dokkum et al. 2004; Forster Schreiber et al. 2004). 

The reddening of the line-emitting gas can be compared to 
that of the stellar continuum emission. Daddi et al. (2004) 
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derive E(B — V) = 0.4 for the continuum of CDFS-695, from 
SED fits. This difference between the reddening of the line- 
emitting gas and the continuum is consistent with previous re- 
sults for nearby star burst galaxies (Calzetti 1997) and DRGs 
(van Dokkum et al. 2004). 



C(B-V)<0.2 0.2<C(B-y)<0.8 E(B-V) > 0.8 




0.1 1 
[Nil] 6583 /Ha 

Fig. 2. — Diagnostic diagram for spectral classification. Small dots 
are nearby galaxies from the SDSS. Galaxies to the left of the dotted 
line are normal star forming galaxies (Kauffmann et al. 2003). The red 
symbol is CDFS-695, which has a Seyfert-like spectrum. The grids 
show shock ionization models from Dopita & Sutherland (1995), with 
and without a precursor region. Histograms show measurements for 
DRGs (van Dokkum et al. 2004) and LBGs at z ~ 2 ([N II] / Ha; Shap- 
ley et al. 2004) and z ~ 3 ([O III] /H/3; Pettini et al. 1998; 2001). The 
blue symbol is the LBG MS 1512-cB58 (Teplitz et al. 2000). The 
top panels show SDSS galaxies binned by reddening: the most dusty 
galaxies often have LINER/Seyfert-like line ratios. 



4. SPECTRAL CLASSIFICATION 

Galaxies whose emission lines arise from photo-ionization 
in H II regions follow a well-defined sequence in diagrams in- 
volving the line ratios [O III] / H/3, [N II] / Ha, [S II] / Ha, and 
[O I] / Ha, largely because of strong correlations between ion- 
ization parameter and metallicity (e.g., Veilleux & Osterbrock 
1987). The diagnostic diagram of [O III] A5007/H/3 versus 
[Nil] A6583/Ha is particularly useful, as these lines are rel- 
atively strong and reddening corrections are small. 

Figure 2 shows the distribution of nearby galaxies in this 
diagram using data from the Sloan Digital Sky Survey DR2 
(Kauffmann et al. 2003; Tremonti et al. 2004). The line ratios 
were corrected for reddening, and only the ~ 124,000 galax- 
ies for which all four lines are detected at > 2a are shown. The 
dotted line shows the separation between H II type emission and 
LINER/Seyfert type emission as adopted by Kauffmann et al. 
(2003). The emission lines of galaxies to the right of this line 



are affected or dominated by other sources of ionization than 
massive stars - in most cases an active nucleus. 

The solid red symbol in Fig. 2 shows the reddening-corrected 
lineratios of CDFS-695. Its [Nil] /Ha ratio of ~ 0.6 and 
[O III] /H/3 ratio of ~ 6 place the galaxy firmly in the Seyfert 
region (see, e.g., Veilleux & Osterbrock 1987). 

TABLE 1 

Emission Lines 
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[Ne III] 
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H/3 
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[O m] 


4959 


43 


10 
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[O m] 
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24 
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[Ol] 
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< 8 
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<0.7 




[Nil] 


6548 


27 


6.4 


2.4 
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Ha 


6563 


99 


24 


9.0 
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[Nil] 


6583 


60 


14 


5.4 


0.1 


[S II] 


6718 


19 


4.5 


1.7 


0.2 


[S II] 


6730 


17 


4.1 


1.6 


0.2 



Flux in units of 10~ 17 ergs s -1 cm~ 2 . 
b Luminosity in units of 10 42 ergs s~ 1 . 
c Relative uncertainty in flux and luminosity. 

5. IONIZATION MECHANISM 

The most straightforward interpretation is that the line emis- 
sion arises from ionization by the power-law continuum of an 
AGN. As the line widths are only w 250kms~ 1 the galaxy 
would be a Type II quasar (QSO) (e.g., Stern et al. 2002), with 
the broad line region hidden by dust. However, other expla- 
nations are also possible. In particular, ionization due to fast 
radiative shocks naturally produces LINER or Seyfert-like line 
ratios (see, e.g., Dopita & Sutherland 1995). 

Whereas it is difficult to unambiguously identify a Type 
II AGN in the rest-frame optical, it can be straightforward 
in X-rays, radio, and/or the rest-frame ultra-violet (UV). The 
rest-UV spectrum is especially important: Seyfert II galaxies 
(and known type II QSOs) typically have featureless continua 
with strong emission lines of Lya and high-ionization species 
such as Nv, Civ, and [Neiv] (Kinney et al. 1991). Daddi 
et al. (2004) obtained a high quality rest-frame UV spectrum 
of CDFS-695. The spectrum is similar to those of normal star 
forming galaxies, and shows no AGN features: Lya is in ab- 
sorption, and familiar rest-frame UV absorption lines such as 
Si IV A 1393, 1402 are strong. The star formation rate implied 
by the de-reddened UV continuum is ~ 500M Q yr _I , in good 
agreement with the star formation rate implied by its (soft) X- 
ray flux and its radio emission (Daddi et al. 2004). In summary, 
there is no evidence for the presence of an AGN other than the 
unusual rest-frame optical line ratios. 

Although we cannot exclude a hidden AGN (or an AGN that 
has recently turned off) as the ionization source, shock ioniza- 
tion is a plausible alternative. The solid line in Fig. 2 shows 
shock models from Dopita & Sutherland (1995). The "shock + 
precursor" models predict similar line ratios as those observed. 
In these models the high [N II] / Ha ratios are produced by the 
shock itself, and the high [Olll] / H/3 ratios arise from photo- 
ionization of preshock regions by extreme UV and X-ray pho- 
tons generated in the shock. 

The shock interpretation is supported by the behavior of the 
emission lines at large distances from the nucleus. As shown in 
Fig. 3 the lines can be traced to ~ lOkpc from the nucleus, and 
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the line ratios remain Seyfert-like: we find [N II] / Ha « 0.3 and 
[O III] /H/3 w 10 in the outer regions of the galaxy. The null or 
inverted gradient in [O III] / H/3 is difficult to explain if photo- 
ionization from the nucleus is the cause. 

Furthermore, there is direct evidence for a strong galactic 
wind which can provide the mechanical energy for the shocks: 
the velocity difference between rest-frame UV absorption lines 
(Daddi et al. 2004) and rest-frame optical emission lines im- 
plies a wind velocity of ~ 300kms~', and the width and tilt 
of the rest-frame optical emission lines demonstrate gas mo- 
tions of ;> 300kms~'. A natural origin of the wind is the 
combined effects of supernovae and winds from massive stars. 
Local analogs to CDFS-695 could be NGC 3079, which has a 
windblown "superbubble" of diameter 1 . 1 kpc displaying sim- 
ilar line ratios as CDFS-695 (Veilleux et al. 1994), and NGC 
1482, in which the entire (bi-polar) wind appears to be shock- 
excited (Veilleux & Rupke 2002). 




Fig. 3. — Two-dimensional spectrum showing position along the y- 
axis and wavelength along the x-axis, in the vicinity of 1.6^im and 
2.1 fjm. Negative images are caused by the dither pattern within the 
6" slit. The emission lines are clearly extended in the spatial direc- 
tion, and show a pronounced velocity gradient. Line ratios remain 
Seyfert-like even at large distances from the center of the galaxy. 



6. DISCUSSION 

Whatever the cause of the Seyfert-like line ratios, the impli- 
cation is that the interpretation of the rest-frame optical emis- 
sion lines in this galaxy is more complicated than is usually as- 
sumed. Metallicity indicators such as [N II] / Ha, the /?23 index, 
and the [O III] / [N II] index (e.g., Pettini & Pagel 2004) would 



produce erroneous results. Using only the //-band spectrum 
we would derive a low metallicity, and using only the Zf-band 
spectrum we would infer a high metallicity. The Ha kinematics 
are probably also influenced by the outflow and may not trace 
the potential (Veilleux et al. 1994). Finally, the Ha luminos- 
ity can be shock-enhanced and cannot simply be converted to 
a star formation rate. The extinction-corrected Ha luminosity, 
Lhu ~ l.Ox 10 44 ergss~', would imply a star formation rate 
of ~ 800M Q yr _1 (Kennicutt 1998). However, the luminosity 
produced in the shock is of the same order as the observed lu- 
minosity: L Hq ~ 2.4 x 10 26 noAV s 3 , with n Q ~ 10 2 cm~ 3 , V s ~ 
300km s" 1 , and A - 10 8 pc 2 (Binette et al. 1985). 

A critical question is whether CDFS-695 is exceptional or 
typical. Its [N II] / Ha ratio is similar to the three other DRGs 
for which these lines have been measured (see Fig. 2), suggest- 
ing LINER/Seyfert like line ratios may be common in red z > 2 
galaxies. In Lyman break galaxies (LBGs) at z ~ 2 this ra- 
tio is usually lower (Erb et al. 2003), even for the subsample 
with K < 20 (Shapley et al. 2004). It is interesting to note that 
the [O III] / H/3 ratios of LBGs at z ~ 3 are very high (Pettini 
et al. 2001; see Fig. 2), although still consistent with the low 
metallicity tail of the distribution of normal galaxies. The only 
LBG for which both line ratios have been measured so far is 
the lensed z — 2.7 galaxy cb58 (Teplitz et al. 2000). It falls on 
the relation defined by normal star forming galaxies in Fig. 2, 
but there is a factor of ~ 5 difference between the N abundance 
inferred from the [N II] emission and from the N I absorption in 
the rest-frame UV (Pettini et al. 2002). This may hint at low- 
level effects of shocks, but other explanations are also possible. 

The effects of shocks and possibly AGN may be more im- 
portant for the dusty DRGs than for LBGs. The top panels 
in Fig. 2 show the distributions of SDSS galaxies binned by 
reddening. Interestingly, a large fraction of the most obscured 
galaxies show LINER and Seyfert-like line ratios, qualitatively 
consistent with earlier results for luminous infrared galaxies 
(Heckman et al. 1987). Submm galaxies at high redshift ex- 
hibit a large range in [N II] / Ha ratio (Swinbank et al. 2005), 
and it will be interesting to see where they fall on diagnostic 
diagrams like Fig. 2. More generally, it is clearly important to 
measure combinations of line ratios for a large, homogeneous 
sample of z ~ 2.5 galaxies, and to investigate correlations with 
dust content, kinematics, stellar mass and nuclear activity. As 
demonstrated in this Letter, Gemini/GNIRS is very well suited 
for such studies. 



We thank Max Pettini for insightful comments on the 
manuscript. 
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